Generating Function Formula of Heat Transfer in Harmonic Networks 
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We consider heat transfer across an arbitrary classical harmonic network connected to two heat 
baths at different temperatures. The network has TV positional degrees of freedom, of which Nl 
are connected to a bath at temperature Tl and Nr are connected to a bath at temperature Tr. 
We derive an exact formula for the cumulant generating function for heat transfer between the 
two baths. The formula is valid even for 7^ Nr and satisfies the Gallavotti-Cohen fluctuation 
symmetry. Since harmonic crystals in three dimensions are known to exhibit different regimes of 
transport such as ballistic, anomalous and diffusive, our result implies validity of the fluctuation 
theorem in all regimes . Our exact formula provides a powerful tool to study other properties of 
nonequilibrium current fluctuations . 
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INTRODUCTION 

Nonequilibrium systems typically generate currents of 
mass or energy. Understanding the general features of 
currents and their fluctuations is one of the main goals 
in nonequilibrium statistical physics. In this context 
the fluctuation theorem (FT) is a remarkable discovery 
[l|, Q ■ This has been theoretically [ll-[l] and experimen- 
tally demonstrated [H'l in many systems. These studies 
have pointed to the significance of the large deviation 
function (LDF) of the current and the related cumulant 
generating function (CGF) in understanding nonequilib- 
rium steady states. Some other interesting related de- 
velopments include the study of Bertini and co-workers 
who introduced a hydrodynamic fluctuation theory 
to study large dynamic fluctuations in steady states and 
that of Bodineau and Derrida 0, Q who conjectured an 
additivity principle for the LDF and CGF of current, 
from which one can predict the quantitative behavior of 
higher order correlations of currents. 

In the context of transport, most analytic results on 
LDFs and FTs are on systems where the bulk dynam- 
ics is stochastic, such as the simple exclusionprocess, 
zero range process, Brownian motors, etc. di, It 
is of general interest to consider and develop these ar- 
guments for systems with bulk Hamiltonian dynamics. 
However, exact analysis here is generally difficult. For 
example even the problem of demonstrating Fourier's law 
of heat conduction in a deterministic system has proved 
to be difficult and has led to many surprises lol- 12|. 
Heat conduction in harmonic lattices is one exception 
where many nonequilibrium properties can be precisely 
discussed. The average heat current, the main focus 
of work so far, is given by a Landauer-like formula in 
terms of phonon transmission coefficients [l^l • Using the 
transport formula it was recently demonstrated [14] nu- 
merically that disordered harmonic crystals in two and 



three dimensions can exhibit different regimes of trans- 
port such as ballistic, localized, anomalous and diffusive. 
Given that this simple deterministic model exhibits var- 
ious regimes of transport, it is of interest to study the 
generic features of nonequilibrium current fluctuations 
in this system. 

In this paper, we derive the general formula of the 
CGF for heat current in an arbitrary harmonic lattice 
connected to two heat baths, which provides a basis to 
explore generic features of current fluctuations. Consider 
heat transfer through a system from a bath at tempera- 
ture Tl to a bath at temperature Tr. Let Q be the heat 
transferred from the left reservoir to the system during 
measurement time r. In general, the distribution of heat 
P{Q) has an asymptotic form P{Q) ^ e'^'''^^^ at large r 
where h{q = Q/t) is the LDF. The CGF ^{\) generates 
cumulants of the heat transfered and is deflned through 
the relation: (e^'^) _ e^'^(^). The LDF h{q) and the 
CGF ii{X) are connected through the Legendre trans- 
form /x(A) = maxg [Xq -\- h{q)] . Properties of heat cur- 
rent fluctuations are contained in h{q) or equivalently in 
/i(A) and various results such as the steady state FT and 
the additivity principle conjecture can be stated in the 
framework of either the LDF or the CGF. For heat con- 
duction, the steady state FT of Gallavotti-Cohen (GC) 
,2] implies the symmetry relation /x(A) = /i(— A — A/3), 
where A/3 — I/Tr — l/T^, and is referred to as the GC 
symmetry. There are examples where the symmetry of 
/i(A) does not imply the FT [l5[ . However, the CGF and 
its symmetry property themselves provide important in- 
formation on current fluctuations. Among these, one of 
the most interesting consequences of the symmetry rela- 
tion is that it leads to the standard linear response results 
such as Onsager reciprocity and Green-Kubo relations 
[3, [3], and in addition makes predictions of responses 
in the far from equilibrium regime [l7l - [l9l |. So far, for 
Hamiltonian systems, the CGF has been analytically ob- 
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FIG. 1: (color online). A schematic picture of the harmonic 
network connected to reservoirs. Positional degrees of free- 
dom is labelled by i = 1, • ■ • , A^. The points at £ = i £ L and 
r = i £ R are attached to the left and right reservoirs. 




and for 
. Here 



tained for one and two particle systems 
a one-dimensional quantum harmonic chain [2 
we obtain a general formula for the CGF of a harmonic 
system in terms of the transmission matrix of a phonon 
mode uj from one reservoir to the other. Remarkably, the 
expression is robust regardless of the complexity of the 
network and the number of particles that are attached to 
reservoirs and always satisfies the GC symmetry. 



HARMONIC NETWORKS AND HEAT 
TRANSFER 

We consider an arbitrary classical harmonic network 
with N positional degrees of freedom labelled i — 
1,2, ■ ■ ■ , N of which Nl are connected to a bath at 
temperature and Nji are connected to a bath at 
temperature Tr (see figure 1). As an example, for a 
three-dimensional cubic crystal consisting of atoms 
with vector displacements, and with two opposite faces 
coupled to heat baths, we would have N = 3N^ and 
Nl — Nn — 3N^ . We denote the positional degrees 
and their corresponding velocities by the column vectors 
X = {xi,X2,--- ,xnY and V = {vi,V2,--- ,vn)'^ re- 
spectively. We consider the following general harmonic 
Hamiltonian for the system: 



H=-V^MV+-X^KX 
2 2 



(1) 



where M = diagjmi, m2, • • • ,tojv} denotes the mass 
matrix and K the force matrix for the system. We model 
the heat baths by white noise Langevin equations with 
each variable coupled to a bath having an independent 
Langevin dynamics. Let L (R) refer to the set of Nl{Nji) 
points connected to the left (right) bath. For discrimi- 
nating these points from the bulk points, we use indices 
£ ^ i £ L and r = i E B. . We define the A^-component 



noise vectors = 77 



(L) 



V 



such that 



Vi 



ry^^'' and 



Vr = Vr are non-zero. Also we define the diagonal ma- 



trices 7 = 7^^' such that 



(i) 



7r 



lee = Tu 



7f and 



7r are nonzero. The equations of motion 



for the system are then given by: 
V 

KX -jV + 7j 



X 
MV 



-KX -j^^'^V- 



V 



(L) 



7 



77 



(R) 



(2) 



The noise terms are assumed to be Gaussian white noise 
with zero mean and correlations given by {r]i(t)r]ii(t')) = 

2jeTLSe,e,6{t -t'), {vAt)Vr'it')) = 2jrTRSr,r'S{t - t'), 

and {r]i{t)rir{t')) = 0, where we have set the Boltzmann 
constant to the value one. The initial state at i = is 
chosen from the steady state distribution, and we mea- 
sure the heat Q flowing from the left reservoir into the 
system between the times t — to t — t . We thus have: 
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(3) 



A solution of the linear equations Eq. ([2|) can be obtained 
by introducing the following discrete Fourier transforms 
and their inverses: 

{X{t),V{t),Tj{t)}^ ^ {X{iu^),Viu;r^),v{iOn)}e-'^-' , 

71— — 00 

{X{u;n),V{uj,,),fi{uj^)} ^- I {X(t),V{t),T^(t)} e^-"* , 



where a;„ = 27rn/T. Plugging these into Eq. ([2]), we get: 

V{uJn) = -iLo.aG+{u,,) [r)(^)(^„)-|-77(«)(w„)] 



-G+(cj„) [ Kl^X + iuJnMAV ] 



(4) 



G+{oj,,) = [-Afw2^K-S(^)(c^„)-S(«Hu;„)] ,(5) 

where 'E'-^'^\oj) = iujj'^^'^'> , AX = X{t) - X{0), and 
AV = V{t) - V{0). The matrix G+ is the Green's func- 
tion connecting bulk variables with reservoir properties. 
The noise correlations in the Fourier space are given by: 



(%(w„)%'(w„-)) = 2Si^i,S 

71,-71' leTL It , 
(r)r(w„)?7r'(w„')) = 2br^r'^n-n'lTTB.lT. 



(6) 



Since the noise strength r](iiJn) ^ 0(l/r^/^) and 
AX, AV ^ 0(1) we see that the second term in Eq. (|4]) 
is ^ order smaller than the first and so can be 

dropped. It can in fact be shown that it contributes 
order 1/r corrections to the CGF [2^. We note that 
V*{u!n) = V{—LOn), V* {^n) = vi—ijJn) ■ The heat trans- 
fered Q [Eq. ^] can be expressed in terms of the Fourier- 
modes with n > as Q = Z^^^o J2i[~'^leVi{uJn)vg{u}n) + 
7y<;(w„){;|(w„) -I- {;£(w„)77^(w„)]. We define G~{uj) 

G+{-uj) = [G+{uj)]* and t'-^'^^lj) = ImjE^^'^H^)} = 
uj^i^-R), On using the solution ^ without the second 
term, i.e., f (w„) = -iw„G+(a;„) [ry^^H'^n) + v'^^H^n)], 
we get the expression of Q as: 
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n=0 



-iUn[flt{iOn){Gj,i, - Gg,^ + flr{uJn){G^g - G^g)] fjiiuJn) 
CO / ~* ( \ 



(7) 



n=0 



I 



where tjl {flu) denotes an Nl (Nr) component column 
vector of noise belonging to £ E L {r E R) sites, while the 



Nl + Nfi dimensional hermitian matrix A is given by 



A = 



~iujn[G+]RL-2LJn[G+T(^'>G-]RL -2c^„[G+r(^)G-]flfl 



(8) 



The subscript L and R in the matrices respectively repre- 
sent the space of £ and r-sites. In Eqs.([7]) and the w„ 
dependence of G and T have been suppressed. In what 
follows, the ujn dependence in variables is omitted unless 
it is necessary. In obtaining the (L, L) element of the 
matrix A, we have used the following Green's function 
identity which easily follows from the definition ([5]): 



G+ G = 2iG" (r(^) + r(-"')G+ 
= 2^G+(^(■^' + r(^))G-. 

CUMULANT GENERATING FUNCTION 



(9) 



We now proceed to the calculation of the CGF. The 



characteristic function Z{X) 



) is obtained by using 



the expression in Eq. ^ and averaging over the Gaus- 
sian noise variables whose correlation matrix is given in 
Eq. We get: 







\ A 2Ti, l-^ J 



fiR 



(10) 



where M denotes the normalization factor of the noise 
distribution. By performing the Gaussian integral, one 
obtains the formal expression of the CGF: 

= -logZ(A) =-iVlogdetS , (11) 



B = 1-A 



Tl Tlr 
Trl —Tr 



n>0 

TlI 
TrI 



(12) 



where 
Tl 

Tr 
Tlr 
Trl 



4[G+r(^)G-r(^)]LL, 
4[G+r(^)G-r(«)]flfl, 
_4[G+r(«)G-r(«)]Lfl + 2i[G-r(«)]ifl, 
-4[G+r(^)G-r(^)]flx - 2z[G+r(^)]flL. 



The matrices Tl and Tr are respectively Nl x Nl and 
Nfj X Nr square matrices, and these can be regarded as 
transmission amplitude of energy with the mode a;„ from 
one reservoir to the other. For Nl — Nr, it is known that 
these matrices appear in the Landauer-like formula for 
average current [13]- As clarified later, even for unequal 
case Nl ^ Nr, both of these are transmission matrices 
and enter in the Landauer-like formula. Although the 
physical meaning of Tlr and Trl are not clear, these 
are closely related to Tl and Tr- The relations can be 
revealed by using the relation ^ iteratively. Through 
tedious but straightforward calculations, one finds the 
following nontrivial relations: 



TlrTr = TlTlr, 
TrlTl = TrTrl, 
TrlTlr = Tr{1 — Tr). 



(13) 
(14) 
(15) 



In order to get simple form of the CGF, we need to sim- 
plify the determinant of B in Eq. (ITT]) . The relations 
P3l) - (|15p play a central role in simplifying the determi- 
nant of B and in deriving the final expression of the CGF. 
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We heuristically introduce the matrix C: 



Clr = ^TrTlr - 
Crr — 1 



Tl 



TlTlr 



1 



1 Tr + TrlTlr 



XTlTr) . 



(16) 



(17) 



The advantage of introducing the matrix C is that the 
product BC has a simple form, and this is useful to 
simplify detB given by detBC/detC = det BC/ det Crr. 
With the relations ([T3l) -(fT5 |) . one finds the following form 
for the product: 



BC 



1 - XTlTl 
—XTlTrl 







1 + ^ ) (1 - TrTlTrMX + A/3)) 



(18) 



r 



and hence: 



det B = det (1 - TrTlTr\{\ + A/3)) 



det(l - XTlTl) 
det(l - XTlTr) 



(19) 



Now by taking the singular value decomposition of the 
matrix [{T'^^^fl^G+{T^"-'^Y/'^]LR it can be shown that 
Tl and Tr have the same set of non-zero eigenvalues. 
Hence det(l — XTlTl) = det(l — XTlTr) and on using 
this in Eq. ()19p we get, in the large r limit: 



Ztt 



X Trlog 



1 - T{uj)TlTrX (A + A/3) 



(20) 



where one can use either Tl and Tr, for the transmission 
matrix T{uj), both of which generate the same values of 
current cumulants. This formula for the CGF is the cen- 
tral result of this paper. One can easily check that the 
GC symmetry: fi{X) — /i(— A— A/3) is satisfied. When the 
system is one dimensional and Nl = Nr = 1 , the formula 
reproduces the classical limit of the quantum version of 
CGF [22]. Interestingly, the formula (|20)) is valid even for 
Nl ^ Nr. In this paper, for simplicity we demonstrated 
the derivation for baths with white Gaussian noise. How- 
ever, the formula is also valid for generalized Langevin 
noise with memory kernel, with appropriate definition of 
the matrices S and T. 



{Q'^)c/t by taking the first and second derivatives with 
respect to A: 



{Q)c {Tl - Tr) 



27r 



da;Tr[r(a;)]. 



2vr Jo 



duTr[T\Lu){TR - Tl)'' + 2T{uj)TlTr]. 



Higher order cumulants are also systematically given. 

The transmission matrix can be obtained either ana- 
lytically or numerically. In case of higher dimensional 
regular lattices, the recursive Green's function method 
[23| can be used to efficiently generate the transmission 
matrix and thus evaluate the CGF. The disordered har- 
monic lattice shows different regimes of transport, such 
as ballistic, anomalous and diffusive transports and hence 
our result implies validity of the fiuctuation theorem in 
all regimes of transport. In addition, the formula (|20|) 
can be a powerful tool to explore generic features of cur- 
rent fluctuation in these different regimes. One of the 
most interesting possible application would be a test of 
the conjecture of the additivity principle Q. 

An open problem is the quantum expression of the 
CGF. As in the one-dimensional case [223 > a two-point ob- 
servation protocol to get distribution of heat is necessary, 
and this seems to be a much more complex calculation 
than the one in this paper. 



DISCUSSION 

We have derived an exact formula for the CGF of a 
general harmonic network ((20|) and shown that it satisfies 
the GC symmetry. The formula is expressed in terms of 
the phonon transmission matrix. The CGF can be used 
to obtain the average current {Q)c/t and current noise 
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